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We apply dimensional reduction to the finite temperature U(l)+Higgs theory 
^ ■ and study the properties of the reduced 3-dimensional theory in the broken 
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cally. The two results for the condensate agree well but a 25% difference is 
observed for the scalar correlator, indicating the need for optimized 2-loop 
perturbative results. 
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1 Introduction 



The U(l)+Higgs theory has been studied extensively in four dimensions 
1^, H, but only to a limited extent in 3d @]. The 3d studies are mainly 
motivated by the fact that the theory then is directly the Ginzburg-Landau 
model of superconductivity. That the 3d theory also arises as an effective 
theory of the finite T 4d U(l)+Higgs theory is a strong additional motivation 
to study it. Although the really physical case is the SU(2)+Higgs theory, it 
is nevertheless of great interest to investigate how changing the group from 
SU(2) to U(l) changes the dynamics. This may also help in understanding 
the phenomena in the full SU(2)xU(l) standard model. The purpose of this 
article is to study the 3d U(l)+Higgs theory, mainly in the broken phase and 
for values of the couplings which are appropriate when the 3d theory is an 
effective theory of the finite T 4d U(l)+Higgs theory 0, |^. 

Perturbation theory in U(l)+Higgs model breaks down at the tree-level 
critical temperature, as the expansion parameter diverges when T — >■ Tc. 
Although the lattice differs slightly from the perturbative one, these tem- 
peratures are so close that perturbation theory is not valid at the lattice 
either. At temperatures below but still higher than all the mass scales 
in the theory, the expansion parameter is small and one expects that per- 
turbation theory works rather well there. The main goal of this article is 
to make statements on the accuracy of perturbation theory at intermediate 
temperatures, where the expansion parameter is not too large and also try to 
estimate the temperature, at which perturbation theory is not valid anymore. 
To this end we compute on the one hand analytically in perturbation theory 
and on the other hand numerically with lattice Monte Carlo techniques the 
scalar condensate and the scalar and vector correlators. A comparison of 
the results gives a quantitative assessment of the accuracy of perturbation 
theory. We shall in a later paper |Q study the phase transition properties of 
the theory. 

The accuracy of perturbation theory depends not only on the number 
of the loops but also on the optimization method used. In this paper we 
compare ordinary ^-expansion with another optimized method, namely the 
CW-method. The latter can also be applied when the broken minimum 
is radiatively generated, the first one requires the existence of the broken 
minimum already in the tree potential. 

First we use the dimensional reduction to obtain an effective model for 4d 
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finite temperature U(l)+Higgs theory. The idea of dimensional reduction is 
to integrate out all the non-zero Matsubara modes which are heavy at high 
temperature |Q |^. The effective 3d theory obtained after reduction provides 
an excellent approximation of the original theory when two basic condi- 
tions are satisfied: the coupling constants are small, so that we can relate 
the parameters of the original 4d and the effective 3d models pertubatively 
and T is larger than all the relevant mass scales in the theory P, ^. Di- 
mensional reduction has one major advantage over direct finite temperature 
4d simulations. In particular the dimensionally reduced 3d theory contains 
less mass scales than the original theory (if we integrate over the temporal 
component of the gauge-field, there are two mass scales less), resulting in 
a huge gain in computer simulations. Because all fermions are integrated 
out at the reduction process, one could also study the theories with chiral 
fermions, which are problematic on the lattice in 4d. 

To be able to compare our non-perturbative Monte Carlo data with per- 
turbative results, we have to determine relevant observables accurately. An 
essential ingredient when comparing the perturbation theory to lattice cal- 
culations are the continuum-lattice mapping formulae, which are exact in 3d 



[IC]. Thus the continuum limit can be carried out under controlled condi- 
tions. 



2 Dimensional reduction 

The Lagrangian of the 4d U(l)+Higgs theory at finite temperature without 
fermions is in standard notation 

£ = \f^.F^, + {D^<py {D^<j)) + m^<j)*<f) + X {<f)*<j)f . (1) 

One could also include fermions as this results in the same final effective 3d 
theory but with different parameter values. 

In the following we shall choose the four dimensional coupling g = 1/3 and 
the 4d tree masses mn = —2rn? = 30, 35, 60 GeV and mw = mH\Jg/2X = 
80.6 GeV. We call the arbitary mass parameter ruH the tree-level Higgs 
mass. The choice of these values of rriH is motivated by the fact that for 
smaller values the transition becomes too strong for dimensional reduction 
to be accurate and for larger values perturbation theory becomes increasingly 
unreliable 
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The dimensional reduction of U(l)+Higgs theory is presented in and 
we just collect the results here. After integrating over the scale vrT the typical 
scale of remaining 3d theory is ~ g^T. There is still an adjoint Higgs field 
Aq present in the theory, but its mass rriAQ ~ gT is large compared with the 
scale of the theory. Hence also Aq can be integrated out perturbatively. The 
action which is obtained after eliminating all heavy degrees of freedom is 



5efr[^i(x),03(x)] = j (fx 



63)' + m20*03 + A3 



^393 



The relations between 3d and 4d parameters are: 

gi = g\f^)T, 



A3 = A(/i)r- 



9l 



(2) 

(3) 
(4) 



where we have 
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2 16tc^ \ 



-/3 + 3A3 



gimp 

47r 



^flog^ + 1 
87r2 V 2mn 2 
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m 



D 



f2m = -^gt + 8\3gl - 8X1 

The constant c = -0.348725. 

The mean field dependent masses are 

= g3(p, mj = mK/ig) + 3A30^ ml = mKfis) + X^cf)'^. 



(5) 

(6) 
(7) 

(8) 



The couplings g^ and A3 are 3d renormalization group invariant. When 
propagators are calculated perturbatively we notice that the mass parameter 
ml contains both linear and logaritmic divergences. The appearance of a 
logarithmic divergence introduces a scale fi^ into our theory. We can choose 
this extra scale to be the scale parameter of dimensional regularization in 3d 



3 



MS scheme. The tree parameters must be transformed with renormahzation 
group equations when /is varies if we want to keep the physical content of 
the theory unchanged. 

All the information of the 4d theory is contained in just three parameters: 
a dimensionful gauge coupling constant g^, which defines the typical scale of 
the system and two dimensionless numbers which we define as 

The parameter x is essentially proportional to the ratio of the squares of 
the scalar and vector masses in the broken phase, while y is related to the 
temperature, y (T — Tc)/Tc, so that at tree level y = describes the phase 
transition point. In terms of 4d parameters they can be expressed as 



1 TTiff ^/Sg 



X = ^ - (10) 



9 3 I Stt 



3^3 , , 1\ mjj 




log — \- c + - 



AnVSg 8n^\ ^ 2g 2 J 2g^T^ 

+^(-4 + 8x-8x^)(log^ + c). (11) 

The 4d tree level Higgs mass values mn = 30, 35 and 60 GeV correspond to 
X = 0.0463, 0.0713 and x = 0.254 according to equation (|TUp. 

3 Perturbative calculation of the properties 
of the system in the broken phase 

3.1 3d effective potential 

The tree potential is obtained from the Lagrangian (^ 
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Figure 1: The 2-loop diagrams contributing to the 3d effective potentiaL 



and the 1-loop part of the effective potential in 3d is 

1 



127r 



2mrj. + + m2 



(12) 



The 2-loop potential 


in 3d contains diagrams in Fig. |T] (xl/(167r^)) 


(a) = 


-g^^^Dyys{mT, rriT, mi), 


(c) = 


-A^0^ [3H{mi, mi, mi) + H{mi, m2, 1712)] , 


(/I) = 


-^g'^Dssv{mi,m2,mT), 


(/2) = 


g^mT{nii + 1112), 


(M) = 


-A {?)m\ + 2mim2 + 3m2^ , 



(13) 



where we have used the results and diagram notations of p[ and 

Dssv{rni,m2,M) = (m^ - 2m^ - 2m^) i/(M, mi, ma) 

[mi — m2) 



M2 



-[i7(M, mi, 1712) - H{0, nil, 1712)] 



DvvsiM, M,m) 



+ — (mi + ma) + (mi — ma)^ — Mmima, (14) 



1 77? 

2i7(M, M, m) - -i/(M, m, 0) + —[H{M, m, 0) 
2 

4 

-i7(M, M, m)] + ^ [i7(m, 0, 0) + H{M, M, m) 



-2H{M,m, 0)] 



m 



m \ 



2M 4M2. 



(15) 



H{mi,m2,m^) = log 



1 



(16) 



mi + m2 + ma 



3.2 Choosing the scale 

As discussed previously the only scale dependent parameter in this theory is 
the mass parameter m'^lfis), as the coupling constants do not run. We also 
know that the effective potential is scale independent to a certain order. On 
the other hand for example the tree as well as the 2-loop effective potentials 
on their own are explicitly /X3 dependent, the scale dependence cancels to 
this order when they are combined. The question is now, how to choose the 
scale so that the higher order loop corrections are small and predictions of 
physical parameters can be made. 

In many theories, including this, the problem involves large logarithms. 
The 2-loop effective potential, which carries the logarithms is only weakly /X3 
dependent. One should choose the scale so that the effect of these logarithmic 
terms is minimized and it is quite natural to achieve this by demanding 
that the 2-loop contribution vanishes totally. This is not exactly sufficient 
because the effective potential is not renormalization group invariant, it is 
not a physical quantity. Only the difference V"(0i) — V^(02) is measurable. 
One should then really look at the derivative dV/d(f), which is /is independent. 

In the effective U(l)+IIiggs theory we have three different mass scales 
(mT, m-i, 777.2) and the determination of proper fi^ is nontrivial. Let's first 
assume that 7777- 3> 777i, 7772 so the Higgs mass is so small that we can neglect 
the nil and 7772 terms. Solving then dV2 = we find that the choice of /is 
which minimizes the 2-loop potential is 



where k has to be determined using some other criterion of minimizing the 
effect of 2-loop terms. Here we demand that the value of (p at the critical 
temperature is the same for the 1-loop and 2-loop potentials in the broken 
phase. This parametrization with constant k works rather well, when the 
Higgs mass is small. When all three mass scales (7777^,7771,7772) are included, 
we generalize {^17\} by choosing 



^3 = nrriT, 



(17) 



/^3(0) 



(18) 
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Figure 2: The renormalization group improved 2-loop potential near the 
phase transition point, with itlh = 35 GeV. 

where K{(f)) is fitted numerically at every point (p so that the derivative of 
the 2-loop contribution vanishes. Pictures of potentials with different Higgs 
masses are presented in Figures ^ and 0. The corresponding k(0):s can be 
seen in Figure Note that the value of and therefore also fis^cp) is not 
well determined for very small values of because perturbation theory is not 
valid there. 




3.3 The ground state energy and the scalar condensate 
in the broken phase 

We will now calculate analytically the ground state energy of the theory and 



the scalar condensate using methods described in p, \iA\. 

In perturbation theory the effective potential is derived as an expansion 
where the expansion parameter is the Planck constant h 

v^(0) = x:^"K(0). (19) 

n=0 
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Figure 3: The 2-loop potential with rnn — 60 GeV. 
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Figure 4: Fitted fi;(0):s with = 35 GeV, T = 148.83 GeV ^ % and 
ruH = 60 GeV, T = 226.77 GeV ^ T^. 
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Although the potential is gauge dependent, its value in the broken minimum, 
which gives the vacuum expectation value I'iT) of is gauge independent 
and determined by 

2^ " A A 



n=0 



0. 



(20) 



The aim is now to solve z^(T) from (pO|). 

We look at the so called "classical" regime, where spontaneous symmetry 
breaking occurs already on the tree-level of the potential. This means that 
there is a (non-trivial) solution Vq{T) = —m\/\^ so that 



dVo{<P) 



0, 



(21) 



't>=MT) 

which gives the leading approximation to Vm 



Vrr. 



H(i^o) 



4A,, 



(22) 



If perturbation theory works well we can use this saddle point value as a 
good start (z/(T) = z/o(T) ^ OiK)) and solve ( PPD perturbatively 



N 
n=0 



(23) 



where we have denoted ip = u'^iT) and ip^ = i^q(T). 

Perturbation theory works only when the expansion parameter is small 
enough. When we look at 1-loop and 2-loop potentials separately we notice 
that typical 1-loop terms have m\/2T{ with them and 2-loop terms include 
5^3777.^/(27?)^. The expansion parameter p is the ratio of these 



9l 



2nmi 



93^3 



(24) 



So the "classical" regime approximation is valid when 



ml < 0, |7772| > ^2^3 



(25) 
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Using ( |I9| ) and (^3]) we have at the minimum 



0. 



(26) 



Here primes mean that we take partial derivatives with respect to 0^ and 
evaluate the expression at the value written in parenthesis. When we also 
expand all Vn-s near (po we get 

^o(<^o) + ^o"(<^o)[^<^l + ^V2] 



All degrees of h must vanish separately, which gives 



0. 



^1 



K."(^o)' 



(27) 



(28) 



From now on we will not write down the y^o-term any more, although all 
potential terms are still evaluated at this saddle point. 

The value of the effective potential at the minimum up to 3 loops is 
calculated by expanding (|1^) with (p3D 



V2 + 



(29) 



which gives the ground state energy in the broken phase. With this, we can 
easily calculate the scalar condensate (0*0(/i3)) in the broken phase 



(0>(/^3)) 



dVrr, 



(30) 



Although the 3d-effective potential is calculated only up to 2-loops, we 
can make some predictions of the 3-loop potential term. 

There is a 3-loop order logaritmic /la-dependence in the 1-loop potential 
Vi because of the /za-dependence of mi(0). Because the 3-loop level effective 
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potential is /is-independent, there must be a term in V3 which cancels this. 
So, V3 must include a term 



3a 



dVi 

'<9/i3 



1 ^"3 , 1 
log + - 

rriT 2 

log + - 

rriT 2 



(31) 



and a /xa-dependence of Vi + V3 is of order h^. When we calculate the ground 
state energy the last term of (^) gives a term proportional to l/m2. At the 
saddle point 1712 and this term diverges. To avoid this, there must be 
another compensating term in V3. Calculation of V2V2 + |^iVi gives us this 
term 

32 V 



3b 



9l 



1/rp 



(47r)3 1712 



(32) 



where 



2A3 

9i 



(33) 



There will be other contributions also, which for pure scalar theory have 
been explicitly computed by []T3[. In SU(2)+Higgs theory this unknown 
information was included under one linear term [Q], which can be derived just 
on dimensional grounds. This term depends on masses, and a dimensionless 
constant parameter (3 



3c 



(47r)3 



p. 



(34) 



In U(l)+Higgs theory this kind of linear term is absent [|TT| so the unknown 
part cannot be found using this parametrization. 

The known part of the 3-loop potential V3 can be written in form 



dVl f2m 

drrii IGtt^ 



log + - 

niT 2. 



— U + 3r^ 
32 V 



9i 



(47r)3 m2 



(35) 



To be able to calculate the scalar condensate we must now find Vmin 
explicitly. First we have to calculate v^o and ipi. Using (^) with (|T^ and 
(PD, we get 



A. 



(36) 



11 



(fi = {ArriT + 3mir + j = ^ ( ^ + ^ ) "^t, 

where masses are evaluated at the saddle point 



and 



nil = rrriT = Y^-2m§(//3), m2 = 0, 



9 ' 



drrii 



1 

mi 



Inserting these into the potentials we get 

m|(^3) 



Vi 



127r 



A3 ' 



rriT: 



1. 



where 



Mr) 

n{r) 



167r2 



2 2 



/i(r){log^ + ^')+gi(r) 



2 4 

Now a lengthy computation shows that 



- (2 - + 
(^2 - + ir^^ log(2 + r) - 2r + ^r^ 

--r + -r log(3r) - - Tr . 



4 \ mT/2m 



log + - 

rriT 2 



+/?disc(r 



(47r)3 

1 4 2 

H ( 4 + 3r^ 1 — 

where the divergent m2 part cancels with V3 and 

4 + 3r3 



Pdisc{r) 



32r4(2 + r) 



(37) 

(38) 
(39) 

(40) 
(41) 
(42) 

(43) 

(44) 



(45) 



- 64 - 32r + 64r2 + 8r^ - 36r^ - 40r^ + lOr^ 
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+13r^ + (48r6 + 24r^) log(3) + (64r^ + 32r^) log(r) 
-(32r^ + 16r^ + 32r^ + 16r^) log(l + r) 

+ (128r2 + 64r3 - 96r^ - 48r^ + 32r^ + 16r^) log(2 + r) 



.(46) 



Using ( p5D and inserting these into (|29| ) gives us the vacuum energy density 
in the broken phase up to 3-loops 
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A, 



127r 



my 



167r2 



2 2 

^3 my 



/i(r)(log^ + ^ 
m^- 2 



(An) 



1 1 

log + - 

rriT 2 



(47r)= 



(47) 



To be able to compare the perturbative results with Monte Carlo data we 
must calculate the analytic form of scalar condensate (0*0). This is easily 
done using (pO]). In terms of previously defined parameters x and y the 
condensate is 



(0*0(/i3)) 



9l 



2x 

1 1 



1 

4:71 



rriT 



9l 



/^3 



/i(r) log h gi(r) 



(47r)3 r^mr 



2m 



^1 



^ 2 \ / /is 1 
r + — ] log 

J \ rriT 2 



(48) 



where tut is 



and 



mr = 5-31 



X 



l/(/w) = y + (-4 + 8x 



x")log^. 



(49) 
(50) 
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T 


106.4 


115.8 


125.4 


126.8 


127.5 


(0*0) /T 


2.395(6) 


1.551(8) 


0.807(4) 


0.705(11) 


0.638(5) 


h-exp. 


2.369(6) 


1.515(9) 


0.707(20) 


0.573(27) 


0.495(37) 


cw 


2.392(2) 


1.550(1) 


0.784(2) 


0.672(2) 


0.614(2) 



Table 1: Measured and computed values of (0*0)/T for ttih = 30 GeV. The 
first row gives the lattice value extrapolated to continuum limit. The last 
two rows show the 2-loop perturbative values using both h expansion and 
CW-optimization. The errors on the perturbative values are obtained by 
varying the scale parameter as indicated in the text. 



3.4 Coleman- Weinberg method 

The direct h expansion discussed above for calculating the scalar condensate 
is not valid at the phase transition point because the expansion parameter p 
diverges. It also has higher order /is-dependence, which is significant near the 
critical temperature. We can mainly avoid these problems using another, so 
called Coleman- Weinberg (CW) method which is much less scale dependent 
and can be used even beyond the critical temperature. With the CW method 
one numerically determines the location of the broken minimum fo{T) of the 
renormalization group improved 2-loop effective potential and then calculates 
the ground state energy of the system. Finally one computes the condensate 
using equation (|30|). Because of the singularities at the classical tree-level 
minimum (V3 ~ l/m2 and m2 0), the ground state energy computed by 
CW-method differs from the result of the straightforward ^-expansion. One 
should then redefine the 2-loop potential by including all leading singularities 
in it (Appendix B in fl^)- This should not be a major problem in this case, 



because the difference is of order and we have the effective potential only 
up to fi^ as there are still unknown parts in the ^^-level. 

The results with both of these methods are presented in section 4. The 
remaining /i3- dependence of the CW-method can be used to measure the 
accuracy of the calculation. We fix the scale /is so that the convergence of 
perturbation theory is maximized near the minimum with a method discussed 
in chapter |3]^. We vary this fixed scale used in the calculation of (0*0(/i3)) 
within limits 0.5yUfix < /^a < 2/ifix. To compare these results with the Monte 
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Figure 5: The 1-loop diagrams contributing to the scalar propagator. 
Carlo data we must also run to fis = using the running 



Table shows that the CW-values of the scalar condensate are almost in- 
dependent of yUs and thus more accurate than the results calculated with h 
expansion. 

3.5 Propagators and particle masses 

The perturbative masses of the system are derived as poles of the correspond- 
ing propagators. Here we calculate the scalar and vector propagator up to 
1-loop level. The propagators are 




log—. 



(51) 



(0i(-A;)0i(A;)) 



1 



(52) 



F + mf-n^(P)' 



{A{-k)Aj{k)) = -— -]^^—^i—— + longitudinal. 




A _ EiEi 



(53) 



At 1-loop level the propagators contain diagrams in Fig. ^ and 

The scalar self-energy is the sum of the following terms (Figure ^ 



V — 



A3 [SAoimD + Ao{ml)\ , 



(54) 
(55) 
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Figure 6: The 1-loop diagrams contributing to the vector propagator. 
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9Bo{k'^; ml, ml) + Bo(k'^; m\, m\ 

-25o(P; ml, 0) [k^ + m^l ^ 
+Bo{k^; 0, 0)k^ - 2Ao(mT)^m|} , 

k'^ + 2k'^(m1 + ml) + (m^ — m2 

2 



(56) 



(57) 



4 {^o(A: 

-5o(A;2;m^,0) 
—AQ{m\)m1 + Ao(m^) 
where we have used the integrals 



fc^ + mi 



2\2 



2 2 
mrp — 7712 



Boi^k'^; ml, m: 



m 



(58) 



(59) 



y (27r)'^[p2 + mf][(p + A;)2 + mi] 
i mi + m2 — 

z(A;2)V2 

The vector self-energy contains the following parts (Figure ^ 
= gl[Ao{ml) + A^{ml 



(60) 



(61) 



k'^ + 2{ml + ml) + 



{m\ — m2 



,2\2 



k^ 
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Aoimj) - Aoiml) 

-5o(P;m2,0) 



2 2 

m| — m2 



-Aoiml) -Aoiml) \, (62) 



- 2ml + Q^l - ""^^^ 



A;2 



(A;2 + mj) 



+Ao(m|) 



1 



A;2 



mj — 



A;2 



(63) 



The mean field dependent masses are evaluated at the minimum of the RG- 
improved 2-loop effective potential. 



4 Monte Carlo simulations 

The validity of perturbative results derived in previous section can be studied 
by performing Monte Carlo simulations of 3d U(l)+Higgs theory. Discretiz- 
ing the continuum action (^ in the usual way one gets 

s = PcT^M^-Uu) 
□ 

-PHT.^e [^*ix)U^ix)(Pix + il)] 

X,fJ, 

+ 0*(x)0(x) + ^ (1 - <p*ix)<Pix))' , (64) 



where Ufj,ix) is a gauge field located on the links of the lattice and 0(x) is 
a complex scalar field located on sites. Un is the product of gauge fields 
around a plaquette. This action has three coupling constants Pg, Ph and Pr, 
giving us a three dimensional parameter space. The continuum limit a — 
is described by the point Pg = oo, Ph = 1/3, = 0. 



4.1 Constant physics curve 

To be able to compare the results from our simulations with the continuum 
theory we have to establish the relation between lattice parameters Pg,Ph 
and Pr and continuum parameters x,y and g^. This is done by a constant 
physics curve, which can be calculated using lattice perturbation theory. Due 
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to the fact that the continuum couphngs are dimensionful, an exact relation 
can be found with a 2-loop calculation [|lOl. The relevant relations are 




1 ' 
4A 



G 



y- (2 + Ax) 



^/3g 

4:71 



1 



-4 + 8, 



ix 



167r2 

X {InGpG + 0.09) + 25.5 + 4.6a; 



(65) 



where E = 3.176. We used x = 0.0463 , corresponding to mn = 30 GeV 
according to equation (0), in all simulations and varied the parameter y, 
which corresponds to temperature according to equation ([Tl|) , between — 1 
and —0.025. At the lower end of this range perturbation theory is expected 
to work rather well, while approaching the critical point yc — it should 
break down. 



4.2 Algorithms used 

By writing the Higgs field as 0(x) = R{x)e^'^^^'' and gauge field as U^{x) = 
Qt9^(x)^ obtain the following expression for action 

S = /3GE(l-COs(^^n)) 

X 

-Ph X! -^(^) cos {-(p{x) + 6^{x) + (p{x + p,)) R{x + fi) 

X,fl 

+ 5:i?(x)2 + /?^E(l-i?(x)2)'. (66) 

X X 

We used an overrelaxed heatbath algorithm for the gauge field. The 
heatbath algorithm used for the gauge field is described in [|1^. This local 
heat bath was then combined with standard overrelaxation, where the angle 
is flipped to the other side of the minimum. 

The form of the action would suggest separation of the radial and angular 
modes of the Higgs fleld, in which case one could by change of integration 
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variables eliminate the angle of the Higgs field as a dynamical variable, thus 
possibly saving some computer time. However, the slow modes of the system 
are associated with the radius of the Higgs field and therefore it is essential 
to update radial part as efficiently as possible. The obvious overrelaxation 
algorithm [jl3l for the radial part of the Higgs field doesn't improve the cor- 



relation time significantly. This is due to a rather poor acceptance ratio. 

The best way to cope 0] with this problem seems to be not to sepa- 
rate the radial and angular update of the Higgs field, but rather update the 
cartesian components. This is done by defining two dimensional vectors 



and 



Rcos(p 
Rsimp 

f cos 9 
f sin 9 



(67) 



(68) 



Now the local action 

5'loc = /-RCOS(0 + 0„) + CiR"^ + C2-R'^, 

can be written as 

^loc = V • r + ciR' + C2R\ (69) 

If we now choose the components of the vector r perpendicular and orthog- 
onal to V, denoted by x and y, we get 

^loc = fx + ci{x^ + y^) + C2{x^ + y^)\ (70) 

which suggests the following overrelaxation algorithm: first we update y — >■ 
—y, which corresponds to overrelaxation of the phase of the Higgs field, and 
then we solve the equation (|7^) to get a new value for x leaving the action 
invariant. 

For an exact overrelaxation algorithm one must leave the total path in- 
tegral invariant, in this case this means that if we change x to x' we should 
have 

eM-SUx)) = ' exp(-5ioc(a:')). (71) 
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Following we now adopt an approximate method and choose x' so that 
the local action 5*100 is invariant and to accept this change with probability 



V 



. .dSiocix') .dSiocix) 
mm( — — /- 



dx' 



dx 



I). 



(72) 



In our simulations the acceptance ratio was between 99.5 - 99.9%. This is 
similar to values reported in three dimensional simulations of SU(2)+Higgs 



theory |n 



The simulations were done using Cray C90 at Center of Scientific Com- 
putation in Helsinki. 



4.3 Condensates 

We concentrated mainly on the scalar condensate, for which the continuum 
and lattice normalizations are related by 

= ^{R'). (73) 

The 2-loop perturbative result calculated in MS scheme is given by equation 
(^HD, which must be related to lattice regularization scheme. One gets 



log(6/3G) + C + Y 



9l 



87r2 



(74) 



where S = 3.176, ( = 0.09 and 6 = 1.94. The scale parameter fi has been 
chosen to be equal to g^. 

We can now compare Monte Carlo data with perturbative calculations 
using equation (|7^). To be able to do this we must extrapolate our results 
both to infinite volume and to zero lattice spacing. To achieve the infinite 
volume limit we made sure that our results did not change with increasing 
lattice size. The actual lattices we used were rather large, typically 60^. 
For the continuum limit extrapolation we performed all our simulations with 
several values of (3g- The leading order of finite lattice spacing effects is linear 
in lattice constant. However, at larger lattice spacings we also observed 
higher order corrections and as a consequence we did either a linear or a 
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Figure 7: The continuum limit of scalar condensate i?^ at y = — 1. Also 
plotted are the perturbative results obtained both with and without CW- 
optimization. The 3-loop results refers to known part of unoptimized h ex- 
pansion. 
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Figure 8: Same as Fig. ^but at y = —0.025. 
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Figure 9: Perturbative results of the scalar condensate compared to Monte 
Carlo data. The straight line is the CW-optimized result, which has been 
substracted from the others. 



quadratic extrapolation to continuum as appropriate. Examples of fits can 
be seen at Figures 

Ordinary perturbation theory quickly breaks down as one approaches the 
critical value as can be seen from Figure where we have plotted the 
Monte Carlo data with 1 and 2-loop straightforward ^-perturbation theory 
results. Also plotted in the Figure ^ are the the known part of the 3-loop 
result and CW- improved result. 

The CW-method is in principle applicable even above the critical temper- 
ature and seems to produce reasonable results at y values up to Uc- However, 
the CW-results are consistently below the Monte Carlo data. This discrep- 
ancy is small (typically less than 4%) but statistically meaningful. The sys- 
tematic error of CW method is not sufficient to explain this effect, as can be 
seen from table |l], so the effect is probably due to higher order corrections. 
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Figure 10: Continuum limit of scalar mass at y = —0.05. 
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Figure 11: Continuum limit of vector mass at y — —0.05. 
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Figure 12: Perturbative results of scalar mass compared to Monte Carlo data. 



4.4 Masses 

The masses of scalar and vector particle at finite temperature were measured 
from the exponential decay of suitably chosen correlators. For scalar particle 
we used the operator 

and for vector 

Im (j)*{x)Ui{x)(j){x + i). 

The quality of our data for masses is not as good as the quality of the 
scalar condensate data, as can be seen from Figures |10| and [Tl|, where we have 
plotted the results of the continuum extrapolation of both scalar and vector 
masses. However, the data is accurate enough to make some conclusions on 
the validity of perturbation theory. 

The tree level calculation and the 1-loop values of scalar mass are plotted 
in Figure |l^ with Monte Carlo data. The Monte Carlo data lies between tree 
and 1-loop values and is not consistent with either of them, implying that 
1-loop calculation is probably not sufficient. 

The Monte Carlo data for the vector mass is consistent with both 1-loop 



and tree values, as can be seen from Figure 13. However the statistical errors 
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Figure 13: Perturbative results of vector mass compared to Monte Carlo 
data. 



are so large that it is not possible to make definite conclusions on the validity 
of the perturbation theory. 

The results obtained here are rather similar to those found in 3d SU(2) + 



Higgs theory [Ij] . There it was found that the perturbative 1-loop calculation 
of vector mass agree extremely well with Monte Carlo data. To be able to 
confirm this result in U(l)+Higgs theory we would need to obtain more 
accurate Monte Carlo data for vector mass. 



5 Conclusions 

We have studied the validity of perturbation theory in the broken phase of the 
U(l)+Higgs theory. Firstly, by measuring the scalar condensate we examined 
the improvement of perturbation theory when increasing the number of loops 
and also the effect of the CW optimization procedure. Our results show 
that the 1-loop calculation is accurate only at y values much lower than the 
critical value Uc- At ?/ = — 1, which is the smallest value of y we examined, 
the error was already 13% and at largest value y = —0.025 the error was 
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73%. Even the 2-loop calculation is not accurate near Uc. The error was only 
1% ai y = —1, but 22% aX y = —0.025. However, optimized perturbation 
theory provides accurate estimates at 2-loop level even up to the critical 
temperature. The largest error was 4% dX y = —0.05. This indicates the 
importance of optimized 2-loop perturbative results. 

The 1-loop perturbative calculation of the scalar mass shows deviations 
from Monte Carlo data. The average difference between perturbative and 
lattice results was 25%. We interpret this as a result of the insufficiency of 
a 1-loop perturbative calculation. This interpretation is supported by the 
fact that the 1-loop values of scalar condensate also show similar deviations. 
For the condensate 2-loop computations are available and their inclusion 
improves the situation considerably. For the scalar mass a 2-loop result is 



unfortunately available only for the pure scalar case [|T5|. It would be very 
interesting to have the 2-loop result for the scalar propagator also for gauge 
theories; then perturbation theory could be optimized and probably a much 
better agreement could be obtained. 

The accuracy of our data for the vector mass has to improve before we 
can make any definite conclusions on the validity of perturbation theory for 
that quantity. One possibility would be to use the smearing or fuzzying 
technique [|^, which has proven to be effective in SU(2)-|-Higgs theory [|l] 



at sufficiently high or low temperatures. 
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